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O Abstract 

o 

Let Fn{x) be the partition polynomial Y!k=\Pk{n)x where pkin) is the number of partitions of n 
►j^ with k parts. We emphasize the computational experiments using degrees up to 70,000 to discover the 

Q asymptotics of these polynomials. Surprisingly, the asymptotics of F„(x) have two scales of orders n and 

and in three different regimes inside the unit disk. Consequently, the zeros converge to network of 
curves inside the unit disk given in terms of the dilogarithm. 

1 Introduction 

u 

Let p{n) denote the number of partitions of a positive integer n; that is, the number of ways of writing n 
additively. Euler identified their generating function: 
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^ By inverting this equation, we find an integral form for p{n) as 

I In their celebrated work, Hardy and Ramanujan in 1917 L6| discovered the asymptotics of p{n) by devel- 

oping a new approach "the circle method" to handle the asymptotics of such contour integrals which have 
r. dense singularities over the unit circumference. A basic form of their asymptotics is 

> 

.W = ^."V/^(1+0(1/VS)) 

We emphasize that Hardy and Ramanujan were guided by an extensive list that appears in their paper of 
values of p{n) up to n = 200. 

The partition polynomial F„(x) is a refinement of the partition numbers: 

n 

Pn{x) = Y^pk{n)x'' 

k=l 

where pk{n) is the number of partitions of « with exactly A: parts. For example, F5 (x) =x^ + 2x^ +2x^ + 
X. We note that the partition numbers are recovered as F„{1) = p{n). Further, pk{n) can be computed via the 
recurrence Pk{n) = Pk-\{n — 1) + Pk{n — k). 
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Figure 1: (a) Partition Polynomial Zeros for Degree 200; (b) for Degree 10,000 



Richard Stanley plotted the zeros of F200 and asked what happens as their degrees go to infinity. See 
Figure 1(a). 

Since the zeros of F„ are symmetric about the real axis, many times we restricted our attention to the 
upper half-plane. 

We single out many intriguing features of this plot: 

1. The zeros are clustering about the unit circle; 

2. There is a sparse of zeros in the left half plane; 

3. There are gaps near the angle 27i/3 and at :t = 1 and x = —\. 

4. There may be two scales of zeros for those either outside or inside the unit circle; 

5. There are 0{n) zeros outside the unit circle and 0{^/n) inside. We use ^/n rather than other powers 
of n by inspiration from the Hardy-Ramanujan asymptotics from p{n). 

6. Empirically we expect the order of convergence of the zeros to be C?( 1 / y/n) lITTI . 

We compare these proposed features with the computed zeros of Fioooo in Figure 1(b). At degree 10,000 
all three gaps now disappear while a second family of zeros inside the unit disk appears. This is especially 
noteworthy at a; = 1 since Fn{\) = p{n). Furthermore, the order 0{^) of zeros is confirmed inside the disk. 

These zeros were found using the MPSolve program described in f3l which is our major tool to obtain 
the zeros of high degree polynomials. Its underlying algorithm is based on simultaneous approximation of 
all the zeros and uses the Aberth iteration. One advantage of MPSolve is that it handles integer coefficients 
with arbitrary large precision. This ability is critical since the coefficients of the partition polynomials have 
hundreds of digits; for example, in Figure 2(b) there is a plot of the number of digits of the 25,000 degree 
partition polynomial. Note that the maximal coefficient 168 digits and that the polynomial is unimodal. For 
comparison in Figure 2(a), we included ^500 whose maximal number of digits is 19. 

In working with such polynomials, a natural question is how to test whether the zeros are relatively 
correct. The zeros obtained have, say, 20 digits, yet the coefficients of the polynomial have hundreds of 
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Figure 2: (a) Digits of the Partition Polynomials for Degree 500; (b) for Degree 25,000 



integer coefficients. So there is no possibility of ever checking directly that the zeros through evaluation. 
On the other hand, the zeros do pass a simple numerical test. We observe that the sum of the zeros equals 
the negative of the coefficient of x"^^ in F„(x) which is always 1. 

The sum of the computed zeros are: at degree 6000, -1.0000000000000000014085; at degree 10,000, 
-0.9999999999999999860563; at degree 15,000, -1.00000000000000002067132 + 2 x IQ-^^i; and at de- 
gree 20,000, -1.00000000000000002323212. For F„, the coefficient of x"-^ is 2 and equals the second 
elementary symmetric function of the zeros. At degree 6000, for the computed zeros this evaluates to 
1.99999999999999999377164 + 5.62 x 10-^6/. 

Furthermore, the zeros that MPSolve produces coincide extremely well with the asymptotics we dis- 
covered so we have much confidence in the computation. In our work there was a complementary interplay 
between developing the asymptotics and numerically determining their limiting behavior and their densities. 

The appealing question of determining the limiting behavior of the zeros is intimately connected with 
the analytic problem of finding asymptotic formulas for these polynomials as their degrees go to infinity. 

2 Introduction to the Zero Attractor 

We now formalize what we mean by the convergence of the zeros. Let Z{Fn) denote the finite set of zeros 
of the polynomial F„. Then the zero attractor of the polynomial sequence {F,,} is the limit of Z{F,i) in the 
Hausdorff metric on the compact subsets of C. In Figure 3(a), there is the plot of the full zero attractor as 
well as a closeup of the upper left-hand quarter of the attractor in Figure 3(b). We will describe the curves 
in these plots in Section |4] 

It may be hard to see but there are actually three curves that make up the zero attractor in the upper 
half plane. Furthermore the common intersection point of these curves we call the triple point T. Its polar 
coordinates are r = 0.9783370882,6 = 2.356797156 and its rectangular coordinates arex = —0.692205581 1 
and 3; = 0.6913717463. 

The zeros converge very slowly to this special point. In Figure 4, there are plots of the zeros near the 
triple point T for degrees 400, 5,000, and 50,000: 

In Table 1, we give the total number of all zeros for F„ inside the unit disk which empirically confirms 
that their order inside the disk is 0{^/n). Its first entry is the count of all zeros inside the unit disk. The 
second entry is the count of all zeros that lie in Q2, the closed quarter unit in the second quadrant excluding 
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Figure 3: (a) Full Zero Attractor; (b) Zero Attractor in Left Half Plane 




Figure 4: Zeros for Degrees 400, 5,000, and 50,000 Near the Triple Point T 



the unit circle. The next entries are the counts for the three families of zeros near the three curves seen in 
Figure 3 that in the second quadrant. These families will be identified in Section [4| for now, we call them 
either Family 1, 2, or 3. Note that Family 1 includes zeros that he along the real axis. 
The predicted number of zeros in Table 1 is given by the empirical formula 

# Zeros inside the unit disk ~ 0.9154v^ (1) 

obtained by fitting a multiple of the power law y/n to the zeros to the first four polynomials at degrees 
5000, 10,000, 15,000, and 20,000. Another estimate for these zeros is given in equation Q. 

3 Potential Theory, Statistical Mechanics, and Limits of Zeros 

If /J is a finite measure with compact support on C, then its logarithmic potential function U^{z) is given as 

u^{z)= I logy^ d^i{t). 

jc \z — t\ 
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Table 1 : Number of Zeros Inside the Unit Disk 



Degree 


if ZiCros 


All ZiCros in (J2 


Family 1 


Family 2 


Family 3 


Prediction 


5UUU 


^ A 

64 


3d 


32 


4 


A 




64.0 


1 AAAA 
lUUUU 




c 1 
J 1 


45 


6 


A 



n 1 T 
91. / 


1 CAAA 

15UUU 




61 


53 


1 




1 12.2 


90000 




71 




Q 
o 




1 9Q 


25000 


146 


79 


69 


9 




144.9 


30000 


160 


87 


76 


10 




158.7 


35000 


172 


93 


81 


11 




171.5 


40000 


184 


99 


86 


12 




183.3 


50000 


204 


109 


95 


13 




204.9 


60000 


226 


121 


104 


15 


2 


224.5 


70000 


242 


129 


112 


16 


1 


242.5 



The connection with polynomial zeros comes from the fact that if ai , • • • , a„ G C and A* = ^ Y!k=\ > then 
C/^(z) = i In |P(z) I, where P{z) = (z - ai) • • • (z - 

The hope is that the logarithmic potentials for a sequence of polynomials will converge to the density 
measure on the zero attractor. Classically, one studies the first limit below while we need the following 
second limit below as well because the number of zeros inside the unit disk is 0{y/n): 

From the Yang-Lee theory of phase transitions in statistical mechanics, the limit of ln[F„(x)]/« gives the 
complex free energy and its zeros indicate the presence of phase transitions. 

We found Alan Sokal's formulation [9J to be very helpful because it allows more general normalizations 
than n. Here is his result: 

Theorem 1 [Sokal] Let D be a domain in C, and let zq G D. Let {gn} be a sequence of analytic functions 
in D, and let {an} be a sequence of positive reals such that are uniformly bounded on the compact 

subsets of D. Suppose there does not exist a neighborhood V of zo '^nd a function v on V that is either 
harmonic or else identically —°° such that 

liminfa„ln|g„(x)| < v(x)limsupa;„ln|g„(;v:)| 

for all z€zV. Then zo lies in liminf Zl(g„); that is, for all n sufficiently large, there exists zeros z* ofg„ with 
lim„^ocZ* =Zo- 

Below we will state what the limits in equation Q are and see that they determine on three regions 
inside the disk. This theorem tells us that the zeros accumulate on the boundaries of these regions. Note: 
these limits are very difficult to compute. 
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4 Explicit Description of Zero Attractor 



Let Li2{x) denote the dilogarithm of x; that is, Li2{x) = ^ 
family of (sub)harmonic functions 

Mx) = 3i\^/U2{x)], fk{x) = l3i 



dt. We introduce a related 



^Ju2{x^ 



which are harmonic in different sectors of the unit disk. These functions play a crucial role in the asymptotics 
of the partition polynomials. Next we introduce the curves Q.f inside the unit disk where fk{x) = fe{x) which 
we will informally call ^Li2-curves. The zero attractor consists of portions of Cu, C13, and C23. 

We can now describe how the plots of the zero attractor are done. Since the \/Li2-curves are level 
curves of the real part of a harmonic function, these curves are naturally solutions to an initial value problem. 
For convenience, write Lk{z) = y^LiiC?)/^ and Lk{z)' is just its usual derivative. Then 

dy ^ '-UlLkix+iyy-Liix + iy)'] 
dx 5[Lk{x + iyy-Li{x+iyy]' 



y{xo)=yo- 



The initial condition needs to be found numerically; typically, by a value on the unit circle where 'Li2{e") 
u{t) + iv{t) is known in closed form: 



u{t) 



n=l 



cosnt 3r-6tn + 2TV- 



12 



v(e)= V '^ = - rin(2sin^)^/^, 0<t<2n. 
„ti Jo 2 



Since the \/Li2-curves are very smooth and even Euler's method gives an useful plot for them. 

In establishing the role of the VLi2 -curves for the zero attractor, though, it is better to treat them as 
boundaries for regions of majorization among the functions fi^ inside the upper unit disk. We define 

!^(1) = {x:/i(x) >f2{x),Mx)},^{2) = {x:f2ix)>Mx),Mx)},^{3) = {x:Mx)>Mx),f2{x)}. (3) 

Then 3^(1) consists of portions of the unit circle, C12, C13, and [0, 1] ; 3^(2) consists of portions of the unit 
circle, Cu, C23, and [—1,0]; 3^(3) consists of portions of the unit circle, C13, and C23. 

The realization that the zeros accumulation on the boundary curves made us examine the majorization 
among /i , /2, and /a on the unit circle. Since these are known in closed form, we found that /i {e") dominates 
on the arc ? G [0,613], fsi^") dominates on f G [613,612], and f2{e") dominates ont £ [612,71] where 613 = 
2.066729664 < 27i/3 = 2.094395103, 612 = 2.2536266 < 37i/4 = 2.356194490, and 623 = 2.361704176 > 
37t/4. The following Theorem together with Sokal's result guarantees that there are no further zero families. 

Theorem 2 (a) On the unit disk, fk{x) < max[fi{x),f2{x),fi{x)], for all k> 4. 

(b) For xe^{l),fi {x) > fk{x), for all k>2. 

(c) Forx G !^(2), /2(x) > fk{x),for allk>3 and k=l. 

(d) For X e !^(3), /a (x) > fk{x), for allk>A and for k=\,2. 

The proof is intricate and lengthy and requires subharmonic function theory and special properties of 
conformal maps. The theorem, though, is easy to verify numerically in special cases. 
We next record the values of the normahzed hmits of In |. 



Theorem 3 Outside the unit disk, lim 



n 



Inside the unit disk, lim 

n— »«> 



and lim 

n— »t» 



Inz, Um 

n— >oo 

ln|^«(z)| 



hi|F„(z) 
n 

0. 



ln|z|. 
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By Sokal's theorem, we now know that the 0{n) contribution to the zero attractor is the unit circle. The 
deUcate asymptotic expansions needed to establish the following theorem are outlined in Section[6] 

Theorem 4 1. On region hm i^S^ = ^Lb(z) and Um = 3iy^lM{zj. 

2. On region ^(2), lim M^^M = i./u^ and lim * = \^\Rm^). 



The above theorem shows that while the limit '"1^"^^) is continuous on the unit disk the derivative of the 
limit fails to be continuous exactly on the boundaries of the three regions. Hence we have determined the 
zero attractor inside the unit disk: 

Theorem 5 The zero attractor for the partition polynomials consists exactly of the unit circle together with 
the boundaries of the three regions ^(1), ^(2), and ^(3) inside the unit disk. 

5 Zero Density 

The zero densities are determined by our following theorem: 

Theorem 6 Density Theorem Let G be a conformal map from a neighborhood of an analytic arc C to a 
neighborhood sector S of an arc of the unit circle such that G{C) = S. Let {r„(x)} be a sequence of analytic 
functions on G ^ (S) such that 

Tn{x) = \ + an{x)G'^" {x) +en{x), 
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Table 2: Properties of the Boundary Curves Cu n ^{l), C13 n 5^(1), C23 n ^{2) 



Curve 






C23n5^(2) 


Length 
Density Mass 

Map Gke 
Circular Arc 
Circular Arc 
Relative Weight 


0.9983742022 
2.464527879 

^VLi2W-VLi2(.v2)/2 

[i,Gi2(r)] 

[0,1.996701527] 
0.8591630301 


0.2884481319 
0.367464849 

^^Li2W-VLi2(A-3)/3 

[1.388229082,1.755693930] 
0.1281025124 


0.02220012557 
0.036529069 

[G23(r),G23(^''*'^^)] 

[1.077010447,1.113539516] 
0.01273445753 



where {c,,} is an unbounded sequence of increasing positive numbers and {a„{x)} is a sequence of an- 
alytic functions that satisfy uniformly on G^^{S) as n ^ 00; |a„(x)| > 5 > 0, = o(l), and e„{x) = 
o{an{x)G'^"{x)). Then for any e > all the zeros ofT„(x) lie in C^for n sufficiently large, where Cg is the 
e-neighborhood of the analytic arc C and 

lim - V{5g(,) : z G ZiT„)nC„} ^ 

where ^ is normalized Lebesgue measure on the unit circle restricted to the circular arc G{C). In particular, 
the zero density measure Vzfor the family isVz = fJoG on the arc C. 

To apply the Theorem to {F„{x)} requires special normalizations that can be found in our paper 0. 

Table 2 gives the essential information about the non-circular boundary curves of the domination regions 
^(1), ^(2), and ^(3). The entry "length" means the arc length of the curve while "density mass" means 
the total measure of the curve with respect to the zero density measure Vz- By the Density Theorem, Vz is 
the pull-back of Lebesgue measure on the unit circle under the conformal mapping Gu where M = 12, 13, 
or 23. The entry "circular arc" means the image of Cu under Gke- Finally by "relative weight", we mean 

WU = Vz(Ch)/[Vz(Ci2) +Vz(Ci3) +Vz(C23)]. 

We emphasize that the number of zeros inside the unit disk is determined, of course, by the density 
measure and not by the length of the curves of the zero attractor. Here is some simple evidence for this. For 
a particular polynomial F,,, the number of its zeros near Cte is about Wk.e out of the total number of zeros 
inside the unit disk in the second quadrant. We confirm this in Table 3. There should be about W13 + W23 
or 14% of zeros near curves C23 and C13. At degree 70,000, there are 242 zeros inside the unit disk. About 
14% of them should lie near these two curves or their reflections in the lower half plane. This count yields 
17 expected zeros which coincides with the computed number. 

In Figure 6, there are plots of zeros along Cu and C13 of degree 50,000. On each plot the portion of the 
curve between adjacent boxes have equal zero density measure. This can be understood graphically from 
counting the number of zeros which are marked by a cross between boxes. The horizontal axes in these plots 
are arclength. For C12, = is the origin; for the other two, 5 = is the triple point T. 

Figure 7 contains the density function for these curves. Since the density function for C12 becomes 
infinite at the origin, we truncated its graph at s = 0.1. Furthermore, we see that the density functions for 
C12 and Ci3 are lightest in a neighborhood of the triple point P while heaviest near the origin for C12 and 
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Figure 6: (a) Zeros of Degree 50,000 Along Curve C12 and (b) Along Curve C13 
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Figure 7: Density Functions as a Function of Arc Length for Curves C12, C13, C23 



the unit circle at e'^'^ for C13. This explains the absence of zeros near the triple point T in Figure 5, even at 
degree 50,000. 

We can also explain why there are so few zeros for the C23-family. Its relative weight is 0.0127. Now 
this makes two zeros near C23 more likely for degree 70,000 over 60,000 yet the opposite occurs. At 70,000, 
the zero that does occur is nearly at the center of the curve. A second zero would violate the near uniform 
density. At 60,000, the two zeros that do occur are near the two endpoints which is consistent with the zero 
density. By equation ([T]) and the relative weight of C23, we do not expect a third zero near this curve until 
degree 190,000. 

In Table 3, the predicted number of zeros near either of the two \/Li2-curves C13 and C23 is found using 
the y^-approximation to the total number of zeros inside the unit disk (see equation ([T])) and the relative 
weights of these two curves. 

Here is another way to get a nearly equivalent prediction of the number of zeros inside the unit disk: 
C • ^/n where 

^ ^ Vz(Ci2)+Vz(Cn)+Vz(C23) ^ 0.9130788466. (4) 

71 

Note that C is remarkably close to the value given by the least squares method in equation ([T]). 
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Table 3: # Zeros for Combined Cn and C23 Families 



Degree 


# Zeros 


Prediction 


5000 


4 


4.5 


10000 


6 


6.5 


15000 


8 


7.9 


20000 


9 


9.1 


25000 


10 


10.2 


30000 


11 


11.2 


35000 


12 


12.1 


40000 


13 


12.9 


50000 


14 


14.5 


60000 


17 


15.8 


70000 


17 


17.1 



6 Asymptotics Inside the Unit Disk 

The purpose of this section is to give a flavor of the asymptotics of the partition polynomials needed to 
confirm their behavior already described. Since the proofs are lengthy, we will be brief. 

The generating function for the partition polynomials is a bivariate version of the generating function 
for the partition numbers which gives an integral version of the polynomials: 

We study F„{x) in each of the three regions %,{k), ^ = 1,2,3 by using the Circle Method. We begin by 
writing Fn{x) as a contour integral over a circle of radius r. 

I r P{x,u) ^ 

Fnix) = ——. (p 

^ ^ 271/ J\u\=r M"+l 

/9" 
P{x, exp[2ni{h/k + iz)]) e-^""'"^!^ g-^Ttm^ # 

where the sum is over rational numbers hjk, (6^ ^, ^'l^ are certain subarcs of the circle that contain the point 
^2%ihlk^ The choice of the radius r and the order of the denominators N of the fractions h/k requires great 
care - this is the initial setup of the Circle Method using Farey fractions. See yj, ||2J, or [7|. Both the radius 
r and the order of the Farey fractions N depend on the degree of the polynomial n and the point x inside the 
unit disk, see 131 . 

The next step is to develop the asymptotics of the integrals 4jt 

P(x, exp[27i/(/j Ik + iz)] ) 

This requires a special approximation to the generating function in a neighborhood of the rational points of 
the unit circle. The proof of this is very lengthy and uses Dirichlet L-f unctions and other techniques from 
analytic number theory. See our paper IH for details. 
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For relatively prime integers h < k < N, we introduce the functions Qh,k{s) that come from expanding 
ln[P(A;,e^''^'''/*^+'^)] and are given in the right plane 3i{s) > 1 as the double series 



y.1 Jlnilmh / k 
Qh.k{s):=l^l^—, {Imr. 

m>\l>\ 



I 



We establish that Qh.k{s) has an analytic continuation to C with a unique simple pole at = 1, with residue 



at = 1 is ^Li2(x*^) and its evaluation at is gft,yt(0) = ^ln(l —x^) + Y.i(k\i) j 



1 



Theorem 7 Let h < k be relatively prime integers. Fix \x\ < 1, then there is a neighborhood ofu = e^^''^/^ 
inside the unit disk where the generating function P(x, u) has the factorization 



P{x,e 



2Ki{h/k+iz) 



, 9^(z)>0, 



where 



'i'kiz) 



Wh.k 



iMx^i 

2nP z 



, ihi{z) 



1 

271/ 



Qh,k{s)^{s) {2%z)-'ds 



2k 



ln(l-/)+ £ 



J 



1 



mi) 



I g-2nilh/k _ I 



Theorem [t] is crucial step that allows applying the Circle Method to the integrals Ih^k- Let /max(-^) 
max[/i(x),/2(x),/3(x)] we state the results: 

Theorem 8 (a) For every k>4, there exists £ G (0, 1) so that 

C fmaxi,^) 1 



\h,k\ < 



/n «3/4 



exp [^Jn{{l-^)fk{x) +/maxW)] 



(b) For A: = 1,2, or 3, on %{k) n [/+; 

h,k 



1 1 



1/2 



exp 2y/n 



\/Li2(^ 



1 + 



uniformly on compact subsets. 



The above Theorem shows that /o,i is the main contribution to on region ^(1), /i,2 the main 

contribution to F„{x) on region ^(2), and /i 3 and /2,3 are the main contributions to F,i{x) on region ^(3). 
In our last Theorem, we state this precisely. For convenience we introduce the notations for the upper 
open unit disk and 4 for 



1 1 



/7I?l3/4 



\/Li2(x) 



1/2 



exp 2^/n 



(5) 



Theorem 9 (a) In ^(1) n U^: Fn{x) = +ok^ {o{h)) uniformly on compact subsets K\. 

(b) In ^(2) n U^: F„(x) = (— l)"e"''-2/2 + ojf2 (o(/2)) uniformly on compact subsets K2. 

(c ) In '^Jyj nU+: F„ (x) = (f,-2it'n/3gWi,3 _^ ^-2Kin2/3^w2.3 ) _^ (/^ ) uniformly on compact subsets K3. 
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We can now describe how to normalize the polynomials F„{x) so we can apply the Density Theorem 
(see Theorem|4]l. 

We can apply Theorem 9 to find the limits ;^ln For example, on region ^(1), 



-^ln|F„(x)| = ^ln|/i(e'*""+o(l))| 



Jn 



\\rvn- ln(27l) + \ ln[Li2(;c)] + 2 V« v^Lii W + 4=1"^"''" +o(l)H 29tv^Li^. 
4 4 \ n 



7 Summary 



We gave an experimental exposition of our work on the zeros of the partition polynomials that are a poly- 
nomial version of partition numbers ||5l. Intensive computational studies are required to understand their 
zeros since the number of zeros inside the unit disk is proportional to the square root of their degree. What 
is especially surprising is that three families of such sparse zeros arise inside the disk. Furthermore, zeros 
for one of these families, the C23 family, do not occur until degree 13,000 and a second zero appears only at 
degree 60,000. 

Our development of the asymptotics for these polynomials completely determines the limiting behavior 
of the zeros both as points and their density. This gives an explanation of the rarity of zeros for C23. 

Among the many polynomial families given by a bivariate infinite product generating function, we 
single out one related to plane partitions. It would be of great interest to find the zero asymptotics for the 
polynomial family {2n(-^)} where they are given by the generating function: 



n=\ k=\ 



1 



We write 2„(x) = Lm=i ^m(")-^'" usual. By fl], Chapter 11, or fTO\, the coefficients qm{n) count the 
number of plane partitions of n whose trace is m; that is, the sum of its diagonal entries is m. Our last figure 
is a plot of the zeros of 2ioo(-^)- 




Figure 8: Plane Partition Polynomial Zeros of Degree 100 
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